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ABSTRACT 

We investigate the use of optical photometric variability to select and identify blazars in large-scale 
time-domain surveys, in part to aid in the identification of blazar counterparts to the ^30% of 7- 
ray sources in the Fermi 2FGL catalog still lacking reliable associations. Using data from the optical 
LINEAR asteroid survey, we characterize the optical variability of blazars by fitting a damped random 
walk model to individual light curves with two main model parameters, the characteristic timescales of 
variability r, and driving amplitudes on short timescales <x. Imposing cuts on minimum r and a allows 
for blazar selection with high efficiency E and completeness C. To test the efficacy of this approach, we 
apply this method to optically variable LINEAR objects that fall within the several-arcminute error 
ellipses of 7-ray sources in the Fermi 2FGL catalog. Despite the extreme stellar contamination at the 
shallow depth of the LINEAR survey, we are able to recover previously-associated optical counterparts 
to Fermi AGN with E > 88% and C = 88% in Fermi 95% confidence error ellipses having semimajor 
axis r < 8'. We find that the suggested radio counterpart to Fermi source 2FGL J1649. 6+5238 has 
optical variability consistent with other 7-ray blazars, and is likely to be the 7-ray source. Our results 
suggest that the variability of the non-thermal jet emission in blazars is stochastic in nature, with 
unique variability properties due to the effects of relativistic beaming. After correcting for beaming, 
we estimate that the characteristic timescale of blazar variability is ^3 years in the rest-frame of 
the jet, in contrast with the ~320 day disk flux timescale observed in quasars. The variability-based 
selection method presented will be useful for blazar identification in time-domain optical surveys, and 
is also a probe of jet physics. 

Subject headings: galaxies: active, BL Lacertae objects: general, quasars: general 



1. INTRODUCTION 

Blazars are a relatively rare sub-class of active galac- 
tic nuclei (AGN) in which a jet is aligned along 
the observer's line of sight, leading to the effects of 
relativistic beaming and unusual associated emission 
(jBlandford fc Reed [l978h . Blazars are among the most 
variable extragalactic objects detected in time-domain 
optical surveys , and have strong emission from radio to 
TeV energies (jUlrich et al. 1997). The central engine 
is believed to be accretion onto a super-massive black 
hole, driving relativistic outflows in a collimatcd jet with 
typical Lorentz factors on the order of T ~ 10 that is 
pointed to within angle of the observer. The term 
'blazars' usually encompasses both BL Lac objects and 
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Flat-Spectrum Radio Quasars (FSRQs), which are be- 
lieved to be jet-ali gned Faranoff-Riley Ty pe I and II 
AGN, respectively ( Urrv fc Padovanilll995l ). In this pa- 
per, we adopt this definition. 

The canonical broadband spectral energy distribution 
(SED) of blazars typically includes several main compo- 
nents: (1) a synchrotron peak likely due to tangled mag- 
netic fields in the jet that may extend from the radio to 
the soft X-ray regime; (2) an inverse- Compton peak in 
the X-ray to GeV regime likely due to scattering of syn- 
chrotron or external photons off of relativistic electrons 
in the jet; (3) a possible inverse-Compton component in 
the soft X-rays due to a hot corona; and (4) occasional 
hints of an underlying accretion disk continuum or host 
galaxy emission. The optical and 7-ray observations we 
use in this paper are expected to be dominated by the 
synchrotron and jet inverse-Compton emission, respec- 
tively. 

The strong high-energy inverse-Compton peak in the 
SEDs of blazars causes them to account for the vast 
majority of bright extragalactic 7-ray emitting sources. 
The Fermi Space Telescope has surveyed the 7-ray sky 
since launch in 2008 , and its Large Area Telescope (LAT, 
lAtwood et aT1l2009f) instrument provides by far the deep- 
est survey to date in the 100 MeV - 100 GeV regime. The 
curre nt 2-year LAT source catalog (2FGL, Nola n et al.l 
2012) includes 1873 total sources, ~44% of which are re- 
liably associated with AGN, and an additional ~14% are 
candidate AGN associations. Of the reliably associated 
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AGN, the overwhelming majority are blazars. Approxi- 
mately half of these Fermi blazars are BL L ac objects, 
and half are FSRQs (|Ackermann et al.|[20lTI ). 

The 2FGL catalog is a significant im provement oyer the 
11-month Fermi LAT source catalog (jAbdo et al.H2010ft 
in number of sources, source detection methods, and 
source associations to known objects. However, ^32% of 
sources in the 2FGL catalog still lack reliable associations 
with any known object, many of whic h may be blazars, 
espec ially at high Galactic latitudes ([Ackermann et ail 
2012). Although the angular resolution of the LAT is 
a dramatic improvement over previous all-sky 7-ray sur- 
vey instruments such as EGRET, typical 95% confidence 
error ellipses of sources in the 2FGL catalog are still on 
the order of several arcminutes in size and can contain 
numerous candidate counterparts at other wavelengths, 
making association and identification of 7-ray sources dif- 
ficult. 

Recently, huge efforts have gone into identifying 
Fermi 7-ray sources by positional coincidences of can- 
didate counter parts from multi-wavelength surveys (e.g. 
IStephen et "all 120101: livlaeda et all 12011ft . observations 
of c ontemporaneous v ariability at different wavelengths 
(e.g. IKara et al.1120121 ). and statistical methods based on 
obser ved 7-ray so urce properties (e.g. Ackermann et al.l 
12012ft . Notably, IMassaro et"aTl (|2011l ) applied a vari- 



ant of the well-known method of selectin g AGN by 
their unique colors in the mid-IR (e.g. ILacv et all 
12004: iStern et all 12005ft to data from the WISE survey 
( Wrig ht et al.H2010ft . and use d their selectio n meth od to 
separate blazars from stars. IMassaro et al.l (|2012f ) then 
further utilized their technique to find WISE-detected 
7-ray blazar candidates in the 2FGL, achieving excellent 
completeness, although the efficiency of their method is 
unclear (see discussion in Section 5.1). 

Aside from their distinct mid-IR colors, blazars are 
also unique in their strong variability at nearly all wave- 
length regimes. This will make blazars stand out in 
the flood of time-domain data from current and future 
large-scale optical i maging surveys such as Pan-STARRS 
dKaiser et al.1 12010ft . the C atalina Real-Time Transient 
Survey (brake et al.l l2009|) . t he Palomar Transient Fac- 
tory (PTF. iRauet alj 12009ft, and the Lar ge Synoptic 
Survey Telescope (LSST. Ilvezic et al.l 12008ft . It is thus 
an auspicious time to explore the possibility of selecting 
and identifying 7-ray emitting blazars en masse in Fermi 
error ellipses by their optical photometric variability. 

Blazars have observed variability on timescales from 
hours to years, and their distinctiveness in time-domain 
surveys results from two effects. Firstly, the relativis- 
tic jet strongly beams non-thermal emission along the 
line of sight towards the observer, boosting the lumi- 
nosity by several orders of magnitude. Secondly, rela- 
tivistic Doppler-boosting shortens the observed timescale 
of variability in comparison to that in the rest-frame of 
the outflow. However, the physical mechanisms caus- 
ing the variability are far less clear. Since blazar emis- 
sion is largely dominated by non-thermal emission from 
the jet, it is possible that internal shocks from overtak- 
ing collisions of fluid shells with different velocities along 
the jet produce time-vari able synchrotron emission (e.g. 
iBottcher fc De rmcr 2010), although other emission com- 
ponents such as the accretion disk may also be time- 
variable. 



Other classes of AGN, including the far more numer- 
ous 'normal' quasars (i.e. non-FSRQ, type 1 AGNs), 
are also highly optically variable, although the predomi- 
nance of particular physical var iability mechanism s may 
differ among AGN subclasses. iSesar et al.l (|2007ft have 
shown that ~90% of quasars in Sloan Digital Sky Sur- 
vey (SDSS, lYork et alj I2000D Stripe 82 are variable at 
the >0.03 intrinsic rms variability level (defined in Eq. 
1), and previous studies using large samples of individ- 
ual quasar light curves have shown that quasar optical 
variability is w ell described by a damped random walk 
(DRW) model (IKellv et al l 120091: IMacLeod et all 
Kozlowski et al.l 120101: iButler fc Blooml 120111: IZu et al.1 



2012ft . IMacLeod et all (|2011l hereafter MA11) further 



utilized the DRW model as a tool to separate stars and 
quasars by optical photometric variability in SDSS Stripe 
82, and showed that quasar selection using this method 
can achieve efficiency E > 80% and completeness C — 
90% by variability alone. This is particularly useful for 
selection of quasars in certain redshift regimes, where 
color selection may fail due to overwhelming contami- 
nation from foreground stars. Furthermore, studies of 
the long-term variability characteristics of AGN may also 
provide a unique perspective on accretion disk and jet 
physics. 

The optical variability of small numbers of individ- 
ual blazars has been studied extensively in the literature 
through intensive long-term monitor ing campaigns (e.g. 
iWebbet al J 1 19881 : iCarini et al.lll992ft . However, no con- 
sistent picture of variability has emerged, at least par- 
tially due to the heterogenous nature of these observa- 
tions and the small sample sizes. Blazar light curves from 
a large-scale flux-limited time-domain survey should be 
more homogenous, and the depth of current surveys can 
yield orders of magnitude more blazar light curves than 
possible through targeted monitoring campaigns. How- 
ever, the rarity of blazars places significant constraints on 
the specifications of any time-domain survey in which a 
large number of blazar light curves can be obtained. De- 
spite the success of the SDSS Stripe 82 for time-domain 
astronomy, the depth of the survey (~22.5 magnitude 
in r band) does not compensate for the small sky cov- 
erage (^290 deg 2 ) when it comes to studying blazars. 
We instead employ the recalibrated LINEAR survey, 
which essentially covers the ^10,000 deg 2 SDSS pho- 
tometric footprin t down to ~17 magnitude in r band 
(|Sesar et al.ll2011l hereafter SE11). IBauer et al.1 (|2009aft 
provided a previous study of blaz ar optical va r iability 
in the Palomar-Quest Survey, and IBauer et al.l (j2009b) 
conducted a subsequent variability-based search for new 
blazars within that survey. However, the Palomar-Quest 
survey provided a median of ~5 epochs of observation for 
each object per filter, and could only be used to study the 
ensemble variability of blazars rather than the individual 
sources. By contrast, the ~200 epochs of observation per 
object in the LINEAR survey allow us to directly model 
the light curves of individual blazars in this study. 

The structure of this paper is as follows: In Section 
2, we describe the LINEAR survey and the construction 
of the Variable LINEAR catalog, from which our blazar 
and comparison normal quasar samples are drawn. In 
Section 3, we describe the damped random walk model 
of variability, and our results of modeling the LINEAR 
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light curves of blazars, normal quasars, and stars. In Sec- 
tion 4, we test the efficiency and completeness of blazar 
selection using a DRW model within Fermi error ellipses. 
In Section 5, we discuss this method in the context of its 
future applications, as well as possible implications of our 
work for AGN physics. We summarize and conclude in 
Section 6. 

2. THE LINEAR SURVEY 
2.1. Photometric Data 

All photometric data used in this paper are from the 
archives of the MIT Lincoln Laboratory (MITLL) Lin- 
coln Near-Earth Asteroid Research (LINEAR) survey, 
spanning the period from December 2002 through March 
2008. A review of the original LINEA R near-Earth 
astero id survey program is presented in Stoke s et al.1 
(2000), and the subsequent photometric recalibration of 
the archived data using SDSS to construct the LINEAR 
photometric database is discussed in SE11. We sum- 
marize only the most salient points here. The LINEAR 
survey program used two 1.01m diameter telescopes at 
the Experimental Test Site within the US Army White 
Sands Missile Range in New Mexico, each equipped with 
a 5 megapixel (2560x1960) back-il luminated CCD de - 
veloped at MITLL and described in iBurke etaLl (j!998[ ). 
The recalibrated survey covers ^10,000 deg 2 of sky, over- 
lapping the SDSS photometric footprint, and contains 
over 5 billion photometric measurements of ~25 million 
objects down to a 5<r depth of r < 18 mag. Due to the 
astrometric goals of the original survey and to increase 
S/N, only a single broad filter was used. The cadence 
range spans minutes to years, with a peak at the main 
15-minute cadence, a gap at 8-hours, and a secondary 
peak around 11 days. The median number of good ob- 
servations per object in the full LINEAR catalog is ^460 
within ±10° of the ecliptic plane and ^200 elsewhere. 

SE11 extracted sources from LINEAR imaging using 
fixed-aperture photometry, a nd recalibrated the astrom- 
etry to the USNO-B catalog (|Barron et al.ll2008ft . To re- 
calibrate the photometry using SDSS, SE11 first matched 
LINEAR sources to SDSS DR7 non-saturated, primary 
objects. Using gri photometry of DR7 objects matched 
to LINEAR, SE11 then modeled and calculated synthetic 
LINEAR magnitudes of matched SDSS objects m-SDSS, 
effectively turning the SDSS imaging catalog into a cat- 
alog of LINEAR photometric standard stars. After a su- 
per flat-field correction on each field using these calibra- 
tion stars, a catalog of 5 billion individual point sources 
with recalibrated LINEAR magnitudes vulinear from 
all epochs of observations was positionally clustered into 
25 million objects, each with at least 15 epochs of ob- 
servation. This comprises the full recalibrated LINEAR 
catalog. Checks on the recalibrated photometry by SE11 
show that the median fulinear — m-SDSS residual per 
field has a distribution of about 0.01 mag wide. This 
should should not be confused with the single-epoch pho- 
tometric uncertainties for individual objects in LINEAR, 
which are generally >0.04 mag (see Figure [5]). Fields 
with hilinear — m-SDSS > 0.1 (usually due to variable 
cloud coverage) are removed. 

2.2. The Bright LINEAR Catalog 
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Figure 1. Normalized histograms of the number of light curve 
epochs available for each source in the Bright LINEAR catalog. 
The distributions are shown for 140 blazars (solid), 978 normal 
quasar (dashed), and all other objects (dotted). 

From the full LINEAR catalog, we select a bright sub- 
set catalog suitable for variability science by imposing 
a cut on the minimum number of good observations in 
each light curve of >30, and a LINEAR recalibrated mag- 
nitude cut of 14 < rriLjNEAR < 17 (where photometric 
errors are < 0.11 mag). At magnitudes >17 and <14, sin- 
gle epoch photometric errors rise rapidly due to photon 
noise and saturation, respectively. These cuts provide us 
with ~4.5 million objects in the Bright LINEAR catalog. 

To create a sample of known blazars and quasars in the 
Bright LINEAR catalog, we positionally match the LIN- 
EAR sources to catalogs of known BL Lac objects and 
FSRQs. To find the known BL Lac objects, we match 
to 1,3 71 BL Lac objects in the iVeron-Cettv fc Veron 
2010) catalog, 501 BL Lac objects in the lPlotkin et at 



,2008) catalog, and 637 radio-loud BL Lac objects in 
the lPlotkin et al.l (120101) cata l og (no ne of the 86 radio- 
quiet objects in lPlotkin et al.l t)2<310( ) matched to a LIN- 
EAR object), all using a 3'.'0 matching radius, resulting 
in 101 matched distinct BL Lac objects. Known quasars 
in the Bright LINEAR catalog w ere identified by match- 
ing to the 105,782 quasars in the lSchneider et al.1 (|2010l ) 
SDSS Data Release 7 catalog, also with a 3'.'0 match- 
ing radius, resulting in 1,020 matched distinct quasars. 
We also separate out the FSRQs by positionally match- 
ing these 1,020 qu asars to the CRATES radio survey 
(jHealev et al.ll2007D of 11,131 flat-spectrum (a > -0.5, 
where cx u a ) radio sources using a 3'.'0 matching ra- 
dius, resulting in 42 distinct FSRQs. The 978 non-flat 
radio spectrum quasars remaining are classified here as 
normal quasars (i.e. unlikely to be strongly dominated 
in the optical by jet emission). 

We note that there is ultimately an overlap of 3 ob- 
jects between the BL Lac object and FSRQ samples due 
to double counting objects which appear i n the BL Lac 
object catalogs, the iSchneider et aLl ()2010l ) quasar cata- 
log, as well as the lHealev et al.1 (|2007f ) catalog, which may 
result from the intermittent appearance of broad emis- 
sion lines in their optical spectra. Since our definition 
of 'blazars' here includes both BL Lac objects and FS- 
RQs, we count these 3 objects as blazars (but not normal 
quasars due to their flat radio spectra). In summary, we 
identified a total of 140 blazars and 978 normal quasars 
in the Bright LINEAR catalog. 

The potential of the LINEAR survey for blazar time- 
domain studies can be seen in Figure [TJ which shows 
normalized histograms of the number of LINEAR epochs 
of observation for each of the 140 objects in the blazar 
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Figure 2. Median intrinsic rms variability <r (solid) and pho- 
tometric uncertainty V , ( m L-TiVBAi?) (dash-dot; Equation 7 from 
SE11) of LINEAR objects in the Bright LINEAR catalog, as a 
function of recalibrated LINEAR magnitude. 

sample, the 987 objects in the normal quasar sample, and 
all other objects (mainly stars) in the Bright LINEAR 
Catalog. The histograms are all generally consistent with 
each other, peak at ^200 epochs, and lack objects with 
<30 good observations due to the previously imposed 
cut. 

Before attempting to model the light curve of every 
object, we would first like to compare the general level of 
variability of blazars, normal quasars, and o ther objects 
in the Bright LINEAR Catalog. We follow iSesar et all 
(2007) and estimate the intrinsic rms variability a for 
each light curve, defined as 

a = [£ 2 - ^{m LINE A R ) 2 ] 1/2 (1) 

for all objects with £ > iP(thlineAr) and a — oth- 
erwise, where £ is the rms scatter of the uilinear for 
all observations of each object, and ip( m LiNEAR) is the 
median photometric error of LINEAR objects as a func- 
tion of magnitude found by SE11 (see their Equation 7). 
Since the photometric errors in LINEAR are relatively 
large, £ mainly reflects photometric noise rather than 
variability for the vast majority of objects. Equation 1 
attempts to remove the effects of the median photometric 
error, and thus a is a better measure of intrinsic variabil- 
ity. Figure [2] shows the median a and iP^linear ) of 
all objects in the Bright LINEAR catalog as a function 
of tci LINE ar magnitude. 

We compare the variability of blazars to normal 
quasars and other LINEAR objects, as it is expected 
that blazars should be systematically more variable. Fig- 
ure[3]shows the integrated distribution of blazars, normal 
quasars, and other LINEAR bright objects in the Bright 
LINEAR catalog as a function of the estimated intrinsic 
rms variability a. We estimate that ~36% of blazars, as 
compared to only ~11% of quasars, are variable with a 
> 0.15 mag in this passband. At this level of variability, 
the estimate of a is not significantly affected by uncer- 
tainties in ip( m LiNEAR.)- Only ~4% of the other bright 
LINEAR objects (non-blazars and non-quasars) are vari- 
able above this a > 0.15 mag level, and they are likely 
a combination of variable stars (e.g. RR Lyrae), eclips- 
ing binaries, underestimates of the photometric errors ip, 
source blending, or other photometric issues. 

2.3. The Variable LINEAR catalog 

In order to characterize the photometric variability, we 
need a clean sample of light curves dominated by intrin- 
sic variability rather than photometric errors. Thus, we 
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Figure 3. Fraction of objects with a intrinsic rms variability 
larger than a, as a function of <r, for blazars (solid), normal quasars 
(dashed), and other objects (dash-dot). 

follow the method of ISesar et al.l (|2007l ) to create a clean 
Variable LINEAR catalog as a subsample of the 4.5 mil- 
lion objects in the Bright LINEAR catalog discussed in 
the previous section. 

As a first variability criterion, we require er > 0.1, 
roughly equivalent to selecting sources with a more than 
two times the measurement noise at 14 < uilinear < 
16, and roughly corresponding to the photometric error 
at niLiNEAR ~ 17. Approximately 8% of objects in the 
Bright LINEAR Catalog pass this initial selection cut, in- 
cluding some non-variable objects at the faint end with 
large a due to large photometric errors. To reduce these 
spurious contaminants, we assume that the photometric 
error distribution follows a Gaussian, and place a cut on 
the x 2 P er degree of freedom calculated with respect to a 
weighted mean magnitude and errors from the photome- 
try. A cut of X 2 /dof > 3 leaves 188,745 objects compris- 
ing the Variable LINEAR catalog we use hereafter, each 
typically having ~200 epochs per light curve. 

Matching this variable LINEAR sample to the vari- 
ous catalogs described in Section 2.2 yields 60 blazars 
and 155 normal quasars, including 1 of the 3 overlapping 
objects discussed previously in Section 2.2 that we in- 
clude in both samples. Although the variable LINEAR 
sample cuts out the majority of the known 140 blazars 
and 978 normal quasars in the Bright LINEAR catalog, 
this is not because the majority of blazars and normal 
quasars are non-variable, but rather because the large 
photometric errors in LINEAR overwhelm the intrinsic 
variability of fainter LINEAR objects. This effect can 
be seen in the much larger fraction of blazars (60 out of 
140) that survive these variability cuts, in comparison to 
normal quasars (155 out of 978). Since blazars tend to 
be systematically more variable than normal quasars (as 
we have shown in Section 2.2), their intrinsic variability 
will dominate over photometric errors for a much larger 
fraction of objects. 

3. THE DAMPED RANDOM WALK MODEL 

MA11 designed and tested optical variability-based se- 
lection of quasars in SDSS Stripe 82 by modeling the 
individual light curves of known quasars as a damped 
random walk. We perform a similar analysis of the LIN- 
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EAR light curves of known blazars and normal quasars, 
and show that blazar light curves lie in a distinct re- 
gion of DRW variability parameter space. The DRW 
model statistically parametrizes variability using three 
parameters: a mean light curve magnitude, a charac- 
teristic damping timescale of variability r, and a driv- 
ing amplitude of short-term stochastic variability a (see 
Kelly et al.ll2009t IKozlowski et al.ll20Tl iMacLeod etHI 
201fl|Zu_eiaL||201J, for further discussion). The fitted 
parameters can be expressed in terms of the structure 
function (SF), defined as the rms magnitude difference 
between all pairs of observations in each individual light 
curve as a function of the time-lag. While this model is 
phenomenological, it is likely that the fitted parameters 
reflect the physical processes that cause the variability. 

Following lkozlowski et al.1 (|2010f ) and MA11, we model 
each individual light curve in our sample as a stochastic 
process described by the exponential covariance matrix 

Sij =(7 2 exp(-|t i -tjl/r) (2) 

for each pair of observations at time ti and tj in the light 
curve. This describes a DRW process with a character- 
istic damping timescale r beyond which the structure 
function will asymptote to a constant value of SFqo = 
(t^/t, and a long-term standard deviation of variability 
a = v^SFqo . The short-term driving amplitude of vari- 
ability a = <j-\/2/t determines the rise of SF(A£) for 
At <C t. We model each individual light curve and cal- 
culate <7 and t, alo ng with their likelih ood distributions, 
us ing the method oflPress et al.1 (|1992t ) , its generalization 
in iRvbicki fc Press! (Il992fl. and the fast c omputational 
implementation in iRvbicki fc Press! (|1995l ). The corre- 
sponding structure function for the model for time-lag 
At is 

SF(At) = SF oc (l - e - |At|/T ) 1/2 . (3) 

In fitting each individual light curve with the DRW 
model, we calculate the likelihood of a DRW solution 
LdrW: as we U as the likelihood of a pure white noise so- 
lution L no i se . Light curves that are better described by 
the DRW model than white noise at a 5-sigma level will 
have AL no i se = ln(L DRW I 'L nolse ) > 5. All our FSRQ 
and BL Lac object light curves in the Variable LINEAR 
catalog fit the DRW model at above this 5-sigma level. 
To remove light curves for which the survey length is 
shorter than r (thus leaving t unconstrained), we also 
calculate the likelihood of a runaway timescale, Lrx,. 
Objects for which r is constrained in t he DRW model 
will h a ve AL^ = \n(LpRw I L oo) > 0.05 (MacLeo d et al.l 
120101) . IMacLeod et all (|2010l ) showed that the DRW is 
an excellent fit to quasar light curves, and used cuts on 
minimum AL no i se and minimum ALao to select quasars 
with high quality light curves. In the higher stellar den- 
sity regions of Stripe 82, MC11 used cuts on maximum 
ex, minimum r, and minimum AL no i se to achieve impres- 
sive completeness C — 93% (defined as percentage of the 
total confirmed quasars in the sample which are selected) 
and efficiency E > 78% (percentage of quasar candidates 
selected which are confirmed quasa r s) in quasar selec- 
tion. Furthermore, IKozlowski et al.l (|2010T l showed that 
the method works even in the high stellar density regions 
of the Magellanic Cloud fields 

We calculate the best-fit DRW variability parameters 
of the individual light curves of the 60 blazars, 155 nor- 



mal quasars, and a random sample of 6000 other objects 
(mainly foreground Galactic stars, representative of the 
other > 180,000 objects) from the Variable LINEAR cat- 
alog. To first understand the underlying distribution of 
DRW variability parameters for these different popula- 
tions, we redshift -correct our parameters as outlined in 
Kellv et al.l (12 009) usin g spectroscopic redshifts from the 
Schneider et al.l (|2010l) catalog for all FSRQs and nor- 
mal quasars. For BL Lac objects, redshifts are much less 
accurate (and sometimes not possible) due to the lack 
of strong spectral features. Since our sample of BL Lac 
objects is drawn from a variety of catalogs, there are 
sometimes discordant redshifts reported; in such cases, 
we preferentially adopt the spectros copically deriv e d red - 
shif ts (or the lowe r limit s) from the lPlotkin et al.1 (2008) 
and Plotkin et al. (2010) catalogs. For 3 BL Lac objects, 
there is no redshift reported in any of the catalogs, and 
so for these 3 objects we adopt the mean redshift found 
for all other LINEAR BL Lac objects, z = 0.28. While 
this assumption is not ideal, very few extragalactic ob- 
jects in the LINEAR survey will be at high redshifts due 
to the shallow optical flux limit, so uncertainties on these 
3 redshifts will not strongly impact our results. 

Imposing the restrictions on Aioo > 0.05 and 
AL no i se > 5 leaves us with 51 blazars and 121 nor- 
mal quasars in our sample. Figure U shows the best- 
fit DRW parameters r against SFqo of blazars, normal 
quasars, and other objects (mostly Galactic stars), along 
with normalized histograms of their distributions. Al- 
though we use separate symbols for BL Lac objects 
and FSRQs in Figure |H we do not separate these two 
blazar sub-populations in our subsequent analysis due 
to the small sample size and unknown biases in our 
blazar sample stemming from possible correlations be- 
tween variability and physical properties of the blazars. 
Future surveys yielding larger samples may be able to 
robustly constrain differences between these two sub- 
populations, and provide further insight into their jet 
proper ties. Figure IH is directl y comparable to the anal- 
ysis of IKozlowski et all (|201Q[ ) and MA11, and confirms 
that the structure functions of normal quasar light curves 
tend to have longer characteristic timescales of variability 
r and slightly larger SFqo than stars. More importantly, 
our analysis also shows that blazars have characteris- 
tic timescales r in between those of normal quasars and 
other objects (mainly variable stars), as well as larger val- 
ues of SFqo. The driving amplitude on short timescales 
tr, is related to these two parameters by the relation 
<x = SFoo/y 7 ?, as shown in Figure |U Figure [5] shows 
the distribution for blazars and normal quasars in <r, 
where normal quasars tend to peak at log (a/mag yr~ 1//2 ) 
~ —0.75, while blazars tend to peak at log a ~ 0.0. A 
two-dimensional Kolmogorov-Smirnoff test between the 
distributions of a for blazars and normal quasars in Fig- 
ure[5]gives a p- value of 2.04 x 10 -5 . For the distributions 
of r in Figure HI the p- value is 0.032. 

While MA11 performed quasar selection using this ap- 
proach by placing cuts on maximum a and minimum 
r, we can also see from Figure [4] that a minimum r 
cut would separate blazars from other objects (mostly 
Galactic stars), while a minimum a cut can help sep- 
arate blazars from quasars. However, it is clear from 
Figure 2] that highly efficient and complete quasar and 
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Figure 4. Distribution of variable LINEAR objects in SFoo (asymptotic value of the structure function on long timescales) and r 
(characteristic damping timescale). BL Lac objects (red stars), FSRQs (red circles), normal quasars (black crosses), and contours enclosing 
90% (blue) and 50% (green) of all other objects are shown. The side panels show the projected normalized histograms of the distributions 
in SFoo (top) and r (right) for blazars (BL Lac objects and FSRQs, red solid), normal quasars (black dashed), and other objects (green 
dotted). Dashed yellow lines of constant a for a = —0.75, 0.0, 0.75 are also shown, with a increasing towards the lower right. 
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Figure 5. Distributions of the driving amplitude on short 
timescales a for blazars (solid), and normal quasars (dashed) in 
the Variable LINEAR catalog. 

blazar selection purely by optical variability in the LIN- 
EAR survey is difficult, as there will be either be heavy 
contamination or a low recovery rate. The main cause 
of this is the bright magnitude range (14 < tiiljnear < 
17) of LINEAR, which overwhelmingly probes Galactic 
stars and very few actual AGN (e.g. there are > 180,000 
sources in the Variable LINEAR catalog, only 214 of 
which are known normal quasars and blazars), and so 
deeper survey data containing a substantial fraction of 
extragalactic sources is necessary. However, currently 
available data from deeper time-domain optical surveys 
such as SDSS Stripe 82 have comparatively small sky 
coverage, and do not yield a large enough sample of the 



much rarer blazars. 

4. BLAZAR SELECTION IN FERMI ERROR 
CIRCLES 

Despite the overwhelming stellar contamination in the 
bright magnitude range of the LINEAR survey, our pro- 
posed blazar selection method is still useful with LIN- 
EAR if we restrict our search to the small regions of 
sky associated with the positions of Fermi sources. This 
greatly reduces the stellar contamination, while simul- 
taneously providing a test of a novel method (based on 
long-term photometric variability characteristics) to as- 
sociate 7-ray sources with counterparts at optical wave- 
lengths. We will show the DRW parameters of known 
Fermi AGN (overwhelmingly blazars) , quantify the com- 
pleteness and efficiency of blazar selection using this 
method, and identify new variability-selected candidate 
Fermi blazar counterparts. 

We first positionally match the Variable LINEAR cat- 
alog described in Section 2.3 to the 1873 sources in the 
2FGL catalog, approximating each error ellipse as a circle 
with radius equal to the semi-major axis of the 95% confi- 
dence ellipse. This initial match yields 480 variable LIN- 
EAR objects in 173 Fermi error circles. These 480 vari- 
able LINEAR optical objects in Fermi error circles are 
then positionally matched using a 3'.'0 matching radius 
to t he 929 associated AGN in the 2nd Fermi AGN cata- 
log ()Ackermann et al.ll201lP ). restricting consideration to 
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Figure 6. Distribution of the number of variable LINEAR objects 
lying in all Fermi error circles of radius r (dashed), and only those 
lying in Fermi error circles of radius r which contain a variable 
LINEAR object matched to an AGN in the 2nd Fermi AGN catalog 
(solid), as a function of r. 



those that have Bayesian association probabilities >0.8 
(this association method is described in the Appendix of 
lAbdo et al.ll2010l) . This yields a sample of 47 known (i.e. 
confidently associated) Fermi 7-ray emitting AGN. We 
calculate DRW parameters for the LINEAR light curves 
of all 480 variable objects, and make nominals cut on 
AL noise > 5 and ALoo > 0.05 as in Section 3. 

Due to the bright magnitude range of the LINEAR sur- 
vey, it is likely that the flux limit of the LINEAR survey 
is too bright to actually detect the optical counterparts 
of the faintest Fermi 7-ray sources, leaving us a popula- 
tion of orphan error circles which contain only LINEAR 
contaminants. Since these fainter sources will also pref- 
erentially have larger error circles, we can remove them 
by considering further cuts on the maximum radii of the 
error circles. In Figure [6l we show the distribution of 
the total number of LINEAR variable objects in Fermi 
error circles as a function of error circle radius r, for all 
173 error circles that contain at least 1 LINEAR variable 
object, as well as the distribution for the 47 error circles 
containing a LINEAR variable object that matched to 
an AGN in the 2nd Fermi AGN catalog. Figure [6] shows 
that the number of variable LINEAR objects in error 
circles with known 7-ray AGN counterparts peaks at 2' 
and drops rapidly; beyond 10', there is only 1 Fermi er- 
ror circle that contains a known 7-ray AGN counterpart. 
However, the distribution of the total number of variable 
LINEAR objects in all error circles as a function of radius 
has a long tail to large r, as there are a small number of 
large error circles (r > 15') with large numbers of vari- 
able LINEAR objects, as expected. This suggests that 
a cut on the error circle radius of ~10' is appropriate, 
but to gauge the impact of such a cut, we will in detail 
consider a range of 8', 10' and 12' cuts on the maximum 
error circle size. A cut at r < 10' removes from con- 
sideration only 1 Fermi error circle containing a known 
7-ray AGN, while a cut at r < 8' removes 4. Scatter in 
the negative relation between the optical flux and 7-ray 
error circle size of Fermi blazars may introduce some bi- 
ases in our resulting sample towards blazars with higher 
than average optical to 7-ray flux ratios when we place 
a cut on the error circle size. However, this bias is ex- 
pected to be small, as only of reliably associated 
AGN in the 2nd Fermi AGN Catalog with optical mag- 
nitudes <17 reside in Fermi error circles >10' in radius. 
We note th at although the Fe r mi A GN association suc- 
cess rate of lAckermann et al.l (|2011l ) as a function of r 



will also affect Figure |H1 this effect is secondary to the 
LINEAR flux limit for the high association probabilities 
considered here. 

Figure [7] (left panel) shows the distribution in r and 
SFqo of LINEAR optical counterparts to Fermi 7-ray 
AGN, as well as all other variable objects in Fermi 2FGL 
error circles of r < 10'. We do not apply cosmologi- 
cal redshift-corrections here for blazar selection. Similar 
to the conclusions drawn from Figure [4] there is clear 
separation between 7-ray emitting AGN and other LIN- 
EAR variables in DRW parameter space. The significant 
normal quasar population seen in Figure [4] is minimized 
when looking only in Fermi error circles. As expected, al- 
most all 7-ray AGN recovered in our sample have blazar- 
like variability. Figure[7] (right panel) shows that variable 
objects in error circles with radius r > 10' largely have 
SFqo and r similar to that of stars, and are thus mostly 
contaminants, rather than new blazar candidates which 
were removed by the cut on r < 10'. 

We can select 7-ray AGN by placing cuts on minimum 
r and a. For different log r and log a, we estimate ef- 
ficiency and completeness of 7-ray AGN selection. The 
results for a 10' cut on the error circles are shown in Fig- 
ure [51 Selection cuts at large minimum r and minimum 
g yield efficiencies of as there are no such objects. Fig- 
ure IH1 shows the maximum completeness achievable by 
tweaking the selection cuts on minimum r and <r, as a 
function of efficiency. This shows how the completeness 
changes as a function of efficiency in 7-ray AGN selec- 
tion, for error circles with a maximum radius of 8', 10', 
and 12'. As expected, the efficiency decreases as the 
maximum size of error circles considered increases, at a 
set completeness. 

By jointly maximizing efficiency and completeness via 
the quantity V E 2 + C 2 from the curves in Figure HI we 
can achieve E > 88% and C = 88% for Fermi error circles 
of r < 8' using cuts on log r > 1.35 and log a > —1.20. 
For error circles of r < 10', E > 76% and C = 86% using 
cuts on log r > 1.53 and log a > —0.90, and for error 
circles of r < 12', E > 70% and C = 86% using cuts 
on log r > 1.53 and log a > —0.90. This demonstration 
verifies that we are able to select (in this case, recover in 
double-blind fashion) 7-ray emitting AGN in Fermi error 
circles with high completeness and efficiency by model- 
ing their optical light curves as a DRW, and imposing 
selection cuts on the variability parameters. 

Our definition for completeness is similar to the con- 
vention of many other recent time-domain AGN stud- 
ies (e.g. MA11). The completeness we calculate here 
is specifically for the selection of variable LINEAR ob - 
jects matched to Fermi AGN in lAckermann etldl ([201 ID 
with Bayesian association probabilities >0.8, and the ef- 
ficiency is calculated assuming these associations are all 
correct. However, this efficiency is a lower bound, as it 
will increase if some of the contaminants we encounter 
in selection are actually blazars. This may occur in our 
analysis if some LINEAR variables that pass all DRW 
blazar selection cuts and did not match to a AGN in the 
2nd Fermi AGN catalog may actually be additional 7- 
ray AGN counterparts, not already recognized as such. 
Indeed, this is likely the case for at least one object, lying 
in the error circle of Fermi source 2FGL J1649. 6+5238. 

In the 2nd Fermi AGN catalog, Fermi source 2FGL 
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Figure 7. Left: Distribution of variable LINEAR objects in Fermi error circles of radius r < 10', in SFoo and r, for reliably associated 
Fermi AGN (red stars), and other objects (black crosses). The underlying contours from Figure [4] enclosing 90% (blue) and 50% (green) 
of all other objects in the Variable LINEAR catalog is shown in the background for reference. Right: Same as the left panel, but for Fermi 
error circles >10' in radius. 




Figure 8. Completeness and efficiency of selection of confirmed Fermi AGN in error circles of radius r < 10', as a function of selection 
cuts on minimum t and a . 



J1649. 6+5238 is listed as associated with the radio source 
87GB 164812.2+524023 at a 0.0 Bayesian probability, 
calculated based on the local density of sources from 
catalogs of likely counterparts. This radio source is po- 
sitionally coincident with a LINEAR object at RA = 
16 ft 49 m 24 s .99, Dec = 52°35'15".05, which has DRW pa- 
rameters calculated from its LINEAR light curve of log 
r = 1.936, log a = -0.026, and log SF^ = 0.942, very 
similar to that of confirmed 7-ray AGN in Figure [7] This 
LINEAR source was not counted as a 7-ray AGN in 
our calculations of efficiency and completeness because 
its Bayesian probability of association was below 0.8 in 



lAckermann et all ([201 ID . Furthermore, this LINEAR 
counterpart also did not match to any known blazar in 
the BL Lac object and quasar catalogs used in Section 
2.2. 

Based on these blazar-like DRW parameters calcu- 
lated from the light curve, we suggest this LINEAR vari- 
able object as the plausible optical counterpart to Fermi 
source 2FGL J1649. 6+523, as well as increased confi- 
dence in associating 87GB 164812.2+524023 as the radio 
counterpart. We note that this Fermi source is also asso- 
ciated with this radio source with a 0.82 probab i lity us - 
ing the log AT— log S method in lAckermann et all (| 2 1 If) , 
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Figure 9. Maximum completeness achievable for a given effi- 
ciency of selection, for confirmed Fermi AGN in error circles of 
radius r < 8' (dashed), r < 10' (solid), and r < 12' (dotted) as a 
function of efficiency. 

based on observed properties of candidate radio counter- 
parts. The origin of this disparity lies in the differences 
in the approaches used to calculate the Bayesian and 
logiV-logS' (radio) association probabilities. We opted 
to use the Bayesian probability >0.8 criterion for known 
7-ray AGN, but the use of the log A^-log S* probability 
also recovers this radio source as a likely 7-ray AGN. It 
is not our intention to compare the two statistical ap- 
proaches, but rather to point out that the discrepancies 
in the results of these two methods may lead to different 
results in analysis. In cases where these two methods lead 
to highly discrepant results, our completely independent 
optical variability-based method can provide valuable ad- 
ditional information. 

Finally, we note that our variability-based approach 
may be further applicable to mass identification of unas- 
sociated 7-ray sources using deeper time-domain optical 
surveys. Data from current surveys such as the PTF, 
with single-epoch depth of r ~ 20.6 mag, 30 epochs of 
observation over 8,200 deg 2 of sky, and ~0.01 mag re- 
peatability can be used to identify unassociated Fermi 
blazars using the method we have presented here. Fu- 
ture surveys such as LSST, with single-epoch depth of 
r ~ 24.5 mag over 20,000 deg 2 of sky and ~ 10 3 epochs 
per source can be used to vastly increase the existing 
sample size. 

5. DISCUSSION 

5.1. Efficiency and Completeness 

In Section 3, we noted that blazar selection using the 
Variable LINEAR catalog across the full field of the 
survey (as contrasted to blazar selection only within 
Fermi error circles) was not efficient and complete. At 
first, this may seem to be in stark contrast to the suc- 
cess of quasar variability selectio n using a similar DRW 
modeling approach (e .g., MAIL iButler fc Blooml [20111 : 
IKozlowski et al.ll2010D . This is primarily a consequence 
of the shallow depth of the LINEAR survey rather than a 
reflection of the ultimate performance expected from our 
approach (especially for fainter blazars). At the depth 
of SDSS Stripe 82 (flux limit of g ~ 20.5), objects with 



low-redshift quasar-like colors account for ^63% of all 
variable objects, with Galactic stars making up the vast 
majo rity of the remaining variable objects (|Sesar et al.l 
120071 ). In contrast, at the shallow depth of the LINEAR 
survey (flux limit of r ~ 17), only ^0.08% of objects 
in the variable LINEAR catalog are quasars (see Section 
2.3). As such, quasar selection by imposing appropriate 
maximum a and minimum r cuts (as done in MA11) on 
the LINEAR quasar sample in Figure [4] will not yield E 
beyond a few percent for any value of C > 50% due to 
scatter in the parameters of the enormously larger pop- 
ulation of variable stars. Blazar selection in LINEAR 
will be similarly inefficient. These issues are exacerbated 
by the large errors on the estimated DRW parameters as 
compared to SDSS, caused by the larger photometric un- 
certainties in LINEAR. In any case, current and future 
time-domain surveys such as PTF and LSST will allevi- 
ate both these issues by going many magnitudes deeper 
while providing ~1% photometry. 

As noted in Section 4, our calculated efficiencies of 
blazar selection are lower limits, since calculating the 
true efficiency requires correct identification of every ob- 
ject as either a blazar or contaminant. There are many 
contaminating variable LINEAR objects that fall in the 
Fermi error circles we considered in our analysis in Sec- 
tion 4 that have blazar-like variability (i.e. have DRW 
parameters similar to blazars and are counted as con- 
taminants when jointly optimizing completeness and ef- 
ficiency for known 7-ray blazar selection) . Some of these 
contaminants are certainly blazars, and may actually be 
the 7-ray source. Securing identification of these blazar- 
like variable objects will require spectroscopic follow-up. 
We list in Table [1] the positions, parent 2FGL error cir- 
cle, and DRW parameters of 12 such new candidate LIN- 
EAR variable objects with blazar-like optical variability 
that fall in Fermi 2FGL error circles with r < 10'. We 
also list any previously-a ssociated counterparts o f each 
parent 2FGL source from lAckermann et al.l ([201 1[ ) with 
Bayesian association probabilities >0.8. The LINEAR 
blazar candidates are not positionally coincident with 
these previously-associated counterparts, but are rather 
additional 7-ray blazar candidates in their respective er- 
ror circles. Since it is possible that some Fermi error 
circles may contain more than one 7-ray emitting ob- 
ject, these LINEAR blazar candidates are still worthy of 
spectroscopic follow-up. 

We further note that the efficiencies we have calculated 
in 7-ray blazar selection in Section 4 are conservative, as 
we have considered all variable LINEAR objects lying 
in all Fermi error circles below a certain size. However, 
~10% of sources in the 2FGL are associated with non- 
AGN objects (e.g. pulsars, supernovae remnants, etc.), 
while ~30% of sources are unassociated. We chose to 
include all these error circles in our analysis (and not 
just those containing known Fermi blazars) in a double- 
blind test, as it is a more faithful demonstration of the 
efficacy of this approach. 

We would like to compare the efficiency and com- 
pleteness of our variability-based Fermi blazar selection 
method to other methods. However, a direct comparison 
is difficult; few studies have attempted to use character- 
istic properties of blazars to systematically select Fermi 
blazars in double-blind fashion as we have done. Further- 
more, the completeness and efficiency of those methods 
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Table 1 

J2000 Positions of new candidate LINEAR blazars in Fermi error circles with r < 10' selected by variability 



RA (degrees) Dec (degrees) Fermi 2FGL Source log(r/days) log(o-/(mag yr 1//2 )) 2nd Fermi AGN Catalog Association 



128.5027 


42.2570 


2FGL J0834.3+4221 


1.89 


-0.20 


OJ 45 


153.3687 


34.4883 


2FGL J1013.6+3434 


1.58 


-0.28 


unassociated 


172.2141 


-5.6568 


2FGL J1129.0-0532 


2.20 


0.02 


unassociated 


189.5802 


-20.0237 


2FGL J1238.1-1953 


2.07 


-0.11 


PMN J1238-1959 


219.1272 


23.2824 


2FGL J1436.9+2319 


1.66 


-0.29 


PKS B1434+235 


220.4843 


43.9407 


2FGL J1442.0+4352 


2.00 


-0.26 


BZB J1442+4348 


232.8992 


57.4483 


2FGL J1531.0+5725 


2.50 


-0.04 


unassociated 


232.9009 


57.3541 


2FGL J1531.0+5725 


2.36 


-0.58 


unassociated 


246.9883 


32.4363 


2FGL J1627.8+3219 


2.30 


-0.34 


unassociated 


250.2259 


11.7870 


2FGL J1641.0+1141 


1.89 


-0.20 


MG1 J164058+1144 


252.3542 


52.5875 


2FGL J1649.6+5238 


1.93 


-0.02 


unassociated 


331.7412 


-0.4681 


2FGL J2206.6-0029 


1.80 


0.27 


PMN J2206-0031 



are almost ne ver jointly optimized (as done in Section 4). 
For example. iMassaro et al.l (|2011l ) used mid-IR color se- 
lection to recover known blazars in the WISE survey, and 
extend ed their approach to find new Fermi blazar candi- 
dates (Massar o et al.ll2012f) . By parametrizing the sim- 
ilarity of the mid-IR col ors of WISE sources to WISE- 
detected Fermi blazars, IMassaro et al.1 (|2012| ) imposed 
a mid-IR color-based selection criterion that is able to 
select known Fermi blazars with 87% completeness over 
the full r egion of sky surveyed by WISE in its first year. 
However, IMassaro et al.l (|2012[ ) did not jointly calculate 
the completeness and efficiency of their approach. Their 
mid-IR color selection method has the advantage of the 
all-sky coverage of the WISE survey. However, mid- 
IR color selection is dependent on the degree to which 
non-thermal jet emission dom inates the mid- I R em ission 
of individual blazars. Indeed, Mas saro et al.l (|2011l ) find 
that FSQRs generally have WISE colors closer to nor- 
mal quasars than to BL Lac objects. This may be due to 
significant thermal mid-IR emission in the SEDs of FS- 
RQs, similar to that in normal quasars (e.g. from dust 
emission, see Plot kin et al.ll201~2h . thus making blazar se- 
le ction more difficult. 

IMassaro et al.l (|2012l ) apply their selection method to 
WISE sources in unidentified Fermi error circles to find 
297 7-ray blazar candidates in 156 out of 313 uniden- 
tified Fermi error circles. However, their estimates of 
the efficiency of their method by systematically offset- 
ting the position of each Fermi error circle to assess the 
number of random associations gave 262 false blazar can- 
didates (using the same mid-IR color selection criteria), 
suggesting that the efficiency of this approach may be 
low. Nevertheless, mid-IR color- and variability-based 
selection of blazars are independent and complementary 
methods, with different selection biases. A combination 
of both approaches, using the full WISE survey as well as 
data from current time-domain optical imaging surveys, 
may be highly efficient and complete. 

5.2. Implications for Blazar Variability 

Our finding that optical blazar variability is well de- 
scribed as a DRW process with distinct variability pa- 
rameters may have interesting potential implications for 
blazar jet physics. The SFqq and r DRW parameters 
should both be affected by relativistic effects. For exam- 
ple, the observed blazar characteristic damping timescale 
Tbiz.obs (after correcting for cosmological redshift) should 
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Figure 10. Top: Distribution of calculated kinematic Doppler 
factors of blazars in the Variable LINEAR catalog, assuming that 
the rest-frame characteristic timescale of variability for each blazar 
is equal to that of LINEAR normal quasars (320 days). Bottom: 
Distribution of calculated rest-frame characteristic timescale of 
variability for blazars in the Variable LINEAR catalog, assuming 
a kinematic Doppler factor 8 = 10. 



be shortened in comparison to the interval Tuz,rest in the 
rest-frame of the emitting material due to the relativis- 
tic motion, with Tu z ,rest = 8Tuz,obai where 5 is the kine- 
matic Doppler factor. While Tbi z ,reat is uncertain and dif- 
ficult to measure, we can place constraints on it through 
comparisons to T qso , the characteristic DRW timescale 
of normal quasars (which are not affected by relativistic 
effects), and independent measurements of 6. 

In FigurelH the distributions of r peaks at Tu z ,obs ~ 80 
days for blazars, and r qso ~ 320 days for normal quasars. 
While these distributions have large tails, their relative 
peaks provide an adequate basis for comparison. If the 
underlying variability in quasars ultimately also drives 
blazar variability, then S ~ T qso /Tbi z . bs ~ 4. Figure [TU] 
shows the distribution of this estimated <5 for blazars 



in the Variable LINEAR catalog, assuming r q 
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days. This is smaller than estimates of the typical blazar 
Doppler factor from radio observations of -F errm-detected 
blazars (8 ~ 10- 30. lSavolainen et al.ll2010f ). and may im- 
ply that the variability of the non-thermal jet emission 
in the jet rest-frame (e.g. due to shocks in the jet) occurs 
on longer timescales than the thermal disk variability of 
normal quasars. 

We can instead adopt the observed typical Doppler fac- 
tors for blazars and calculate Tu z ,rest for each individual 
blazar. Ideally, we would use a measured value of 8 for 
each individual blazar in this calculation, but direct mea- 
surements of 8 from the literature are sparse and inho- 
mogeneous. In Figure [TU] (bottom panel), we show the 
distribution of Tuz.rest assuming a reasonable 8 = 10 for 
all blazars. The peak in Figure [10] occurs at Tuz.rest ~ 3 
years, which is longer than the T qso ~ 320 days we mea- 
sure for LINEAR normal quasars. The discrepancy be- 
tween 8 for Fermi blazars implied by radio observations 
and 8 calculated assuming intrinsic variability similar to 
normal quasars suggests different stochastic mechanisms 
driving the variability in the disk and the jet. While a de- 
tailed physical interpretation of this is currently unclear, 
our long-term optical variability results may provide ad- 
ditional constraints on models of blazar jets. 

5.3. Short Timescale Variability of Quasars and 
Blazars 

Observations of intraday variability (also referred to as 
'microvariabilty') of normal quasars have shown a puz- 
zling diversity of properties. Not all quasars exhibit this 
phenomenon, but some show varia bility on the order of 
~0.01 mag on ~1 day ti mescales (iGopal -Krishn a et al.1 
2001 IStalin et all 12005 . [2001 : IGuota fc .Toshil 1200.1 
Mushotzkv et all 20111) . The cause of quasar variabil- 
ity on such short timescales is an unresolved problem, 
although the presence of strong intraday variability is 
correlated with radio-loudness, and it has long been sus- 
pected that a weak jet component may be both the source 
of ra dio emission and th e cause of the stronger variabil- 
ity. iKellv et al.l (|2009D showed that the DRW model 
successfully predicts amplitudes of variability of <0.02 
mag over ^8 hours for quasars, consistent with obser- 
vations from quasar monitoring. This successful predic- 
tion of the short-timescale variability is noteworthy since 
the sampling inter vals of the quasar light curves used in 
IKellv et all (|2009t ) are all >2 days, and thus do not ac- 
tually sample such intraday timescales. 

From the log a distributions for normal quasars and 
blazars in Figure [5] we can calculate the standard devi- 
ation of the expected variability on ~1 day timescales, 
approximated as ay/ At. For quasars, the peak of the 
distribution occurs at ^0.2 mag yr -1 / 2 , thus predicting 
intraday variability ^0.01 mag on 1 day timescales. For 
blazars, the peak occurs at ~1 mag yr -1 / 2 , thus pre- 
dicting variability of ^0.05 mag on 1 day timescales in 
the observed frame. Previous studies of blazar microvari- 
ability have often focused on the most variable blazars, 
and have reported night-to-night variations as h igh as 1.0 
mag (e.g. lCarini eTatl[T99lt [Ghosh et al.ll2000T ). The in- 
traday variability we have calculated for LINEAR blazars 
is from an untargeted (i.e. less biased) survey, and is 
likely more representative of the intraday variability of 
blazars as a whole. 



Despite the large photometric uncertainties in the LIN- 
EAR survey, there are a significant number of normal 
quasars and blazars in Figure [5] which do have extreme 
levels of short-timescale variability, above the photomet- 
ric uncertainty. Nearly all LINEAR objects have photo- 
metric uncertainties below 0.1 mag (and reaching as low 
as <0.04) at <17 mag (see Figured]). Conservatively as- 
suming that variability on 1 day timescales needs to be 
>0.1 mag in order to be detectable above the photomet- 
ric uncertainty implies a log a > 0.28. In Figure [5] 18 of 
119 normal quasars and 10 of 51 blazars in our sample 
are variable above this criterion (and both these frac- 
tions would increase if we more carefully considered the 
photometric uncertainty limit of each individual object). 
Thus for these objects, the intraday variability is above 
the LINEAR photometric uncertainty and the >0.1 mag 
intraday variability is significantly above most previous 
observations of normal quasars. The source of this ex- 
treme variability for a subset of quasars is unknown, and 
warrants additional investigation. 

6. CONCLUSIONS 

We have used light curves from the LINEAR optical 
imaging survey to study the variability of blazars, and to 
compare their variability to that of normal quasars. In 
contrast to many earlier studies of blazar variability, this 
study is the first to use a reasonably large, homogeneous 
sample of individual light curves from a wide-field survey, 
where the distribution of various variability properties for 
blazars and quasars can be calculated and compared. 

In a time-domain survey, blazars are among the most 
violently variable objects detected, and are generally 
more variable than normal quasars. We estimate that 
^36% of blazars but only ~11% of quasars are variable 
at the extreme >0.15 mag intrinsic rms variability level. 
We show that, like quasars, blazar light curves are well-fit 
by the DRW model for variability, but the blazars lie in 
distinct variability parameter-space with higher r than 
stars and higher a than normal quasars. This suggests 
that blazars can be selected with high efficiency E and 
completeness C by imposing minimum selection cuts on 
r and a. 

Due to the overwhelming numbers of variable stars in 
the bright magnitude regime probed by LINEAR, it is 
difficult to select blazars (and quasars) at high efficiency 
and completeness in the full LINEAR sample using op- 
tical variability methods alone. We instead examine a 
more focused test to select 7-ray emitting blazars in 
Fermi error ellipses from the 2FGL catalog. We cal- 
culate DRW parameters for light curves of all LINEAR 
variable objects objects within Fermi error ellipses and 
place minimum selection cuts on r and a. Using this 
approach, we are able to recover the corresponding 7-ray 
emitting AGN counterparts in the 2nd Fermi AGN cat- 
alog with E > 88% and C = 88% for 95% confidence 
error ellipses with semimajor axis r < 8', and E > 70% 
and C = 86% for r < 12'. 

Our 7-ray blazar selection has uncovered a variable 
LINEAR object, coincident with radio source 87GB 
164812.2+524023, in the error ellipse of Fermi source 
2FGL J1649.6+5238. This object was not associated 
with the Fermi source in the 2nd Fermi AGN catalog 
using Bayesian probabilities. Our analysis shows this ob- 
ject has optical variability properties consistent with 7- 
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ray emitting blazars and is likely to be the 7-ray source. 
We find a total of 12 objects with LINEAR variability 
parameters similar to blazars lying in Fermi error cir- 
cles with r < 10' (see Table 1). Confirmation of these 
variability-selected blazar candidates will require spec- 
troscopic follow-up. 

Our results suggest that the variability of the non- 
thermal jet emission in blazars is stochastic in nature, 
with higher amplitudes and shorter observed timescales 
in comparison to the variability of normal quasars, likely 
due to the effects of relativistic beaming. Assuming a 
reasonable Doppler factor of 10, we estimate that blazars 
are characteristically variable on timescales of ~3 years 
in the rest-frame of the emitting material, longer than the 
^320 days characteristic disk flux timescale for quasars. 
This suggests that different physical mechanisms domi- 
nate the observed variability from blazars and quasars, 
likely connected to the jet and to the disk, respectively. 
The fitted parameters imply that blazars have a typical 
intraday variability amplitude of ~0.05 mag, compared 
to ^0.01 mag for normal quasars. Furthermore, there is 
a significant fraction (~15%) of normal quasars which ex- 
hibit large intraday variability of >0.1 mag, detectable 
above the photometric uncertainty. The source of the 
extreme variability of these quasars is unclear. 

We argue that variability-based blazar selection is 
likely to be highly efficient and complete in deeper optical 
time-domain imaging surveys, and that the variability- 
based blazar selection method presented in this paper 
is capable of greatly increasing the number of known 
blazars. Our approach, based on long-term photometric 
variability characteristics of a large sample of individual 
blazar light curves, may also bring a new perspective on 
accretion disk and jet physics. Finally, this work typifies 
the fruitful ancillary science made possible by combining 
data from surveys as diverse in all respects as LINEAR, 
SDSS, and Fermi. 
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